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ABSTRACT 

Newtonian  impact  theory  is  used  to  examine  the  variation  in 
ballistic  coefficient  of  a  coated  cylinder  as  the  thickness  of  the  coat¬ 
ing  decreases. 

It  is  shown  that  when  the  density  of  the  coating  is  less  than 
one -half  that  of  the  core -material,  an  initial  value  of  coating- 
thickness  can  be  chosen  which  results  in  a  steady  increase  of  the 
ballistic  coefficient  as  the  thickness  of  the  coating  tends  to  zero. 

The  largest  such  value  of  coating-thickness  ,  the  correspond¬ 
ing  minimum  value  for  the  ballistic  coefficient ,  and  the  rate  at  w'nich 
the  ballistic  coefficient  rises  to  that  of  the  bare  core  are  shown  to  be 
simple  functions  of  the  particular  density-ratio. 
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INTRODUCTION 


In  order  to  predict  the  atmospheric  deceleration  of  a  cylinder  whose 
axis  lies  in  a  plane  normal  to  the  direction  of  motion  ,  it  is  necessary  to  consider 
the  ballistic  coefficient  ^  =  W/C^^A  ;  W  here  is  the  weight  of  the  cylinder 
and  the  effective  "drag  area"  upon  which  the  dynamic  pressure  acts 

to  oppose  the  motion  of  the  body.  ^ 

Since  relatively  low  ballistic  coefficients  are  associated  with  lightweight 
cylinders  at  large  angles  of  attack,  maximum  deceleration  and  heating  may  occur 
rather  early  following  re-entry.  ^  Survival  of  such  atmospheric  heating  may  be 
enhanced  not  only  by  an  appropriate  choice  of  materials  and  configuration ,  but 
also  by  the  use  of  ablative  coatings. 

We  are  concerned  in  the  discussion  which  follows  only  with  the  variation 
in  ballistic  coefficient  which  occurs  with  mass -loss  when  such  a  protective  coat¬ 
ing  is  used. 


t  CXher  orientations  may  also  be  of  interest ,  and  for  such  cases  /3  should  be 
replaced  in  the  sequel  by  /3  sin^  ^  ,  where  a  is  the  "angle  of  attack" 
measured  between  the  cylinder  axis  and  the  direction  of  flight.  (See  Appendix. ) 
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THE  BALLISTIC  COEFFICIENT- 


If  a  cylindrical  core  of  length  L  ,  radius  r  ,  .  and  density  p  is 
covered  by  an  ablative  coating  of  thickness  R-r  and  density  p'  <  p  ,  then  the 
Newtonian  ballistic  coefficient  of  the  coated  cylinder  when  oriented  normally  to 
the  flight  path  (with  ci  -  ir/2  )  is  given  approximately  by  the  expression^ 


where 


and 


W  2 

P  =CpA  "  l*2Cp’R  +  r  (p-p*)/R] 


W  =  irr\p +ir(R^-r^Lp’ 


CpA  =  (8/3)RL 


(  1) 


Two  limiting  cases  are  of  immediate  interest  for  purposes  of  comparison; 

( 1 )  a  bare  cylinder  of  density  p  and  radius  R  ,  since  it  occupies  the  same 
volume  as  the  coated  cylinder;  and  (2)  a  bare  cylinder  of  density  p  and  radius  r  , 
since  this  is  the  core  which  the  coating  is  intended  to  protect. 

The  first  of  these  has  the  ballistic  coefficient 

J  =  1.2Rp  (2) 

and  the  ballistic  coefficient  for  the  second  is 

P*  =  1.2rp  .  (3) 

t  This  formula  for  the  ballistic  coefficient  and  the  ones  following  are  derived  in  the 
Appendix. 
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where  t<R  (and  ci  =  has  been  assumed  in  each  instance). 


A.  Case  1 

It  is  apparent  from  Formula  (1)  for  the  coated  cylinder  that 

/3-4/3  as  p’-*p  and  as  r-.  R 
Indeed ,  the  defect  from  the  limiting  value  is 

=  1.2  (p_p-)(R2_r2)/R 


=  /3(l-p7p)[l-(r/R)^l 


(4) 


and  under  the  assumption  that  the  coating  is  so  thin  that  (R  +  r  )~2R  ,  .  we  obtain 
the  simple  approximation  that  the  defect  varies  directly  with  the  thickness  of  the 
coating  and  the  difference  in  density  of  the  materials; 


^~2.4(p-o')(R-r) 

or  (5) 

^-^~20(l-pyp)(l-r/R)  . 

A  direct  com^rison  of  the  coated  cylinder  with  outside  diameter  2R  and 
the  bare  one  of  the  same  dimensions  may  conveniently  be  expressed  in  terms  of  a 
"weighting  function" 

g(p’,R)  =  m  (6) 
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where  r  and  p  are  assumed  fixed  . 

Formula  (1)  may  be  rewritten  to  make  the  form  of  g  apparent: 

/3  =  [(p7p)  +  (r/R)^l-p7p)]^  .  (7) 

The  weighting-function  g  is  plotted  in  Fig.  1  for  several  different  values 
of  (p7p)  and  appears  as  an  increasing  function  of  the  parameter  (r/R)  . 

The  penalty  in  ballistic  coefficient  which  a  given  density  and  thickness  of 
coating  may  impose  in  relation  to  a  given  choice  of  core  may  thus  be  assessed 
directly.  (The  densities  of  some  representative  materials  are  given  in  Table  1.) 

However ,  in  viewing  Fig.  1  with  the  dynamic  situation  in  mind ,  it  must  be 
remembered  that  while  g-*l  as  (r/R)-.l  for  each  value  of  (p7p)  >  P 
is  decreasing  as  R-*r  and  consequently ,  so  is  p  . 

B.  Case  2 

Survival  problems  which  arise  from  atmospheric  heating  may  make  the  use 
of  a  protective  coating  attractive  in  spite  of  the  lowered  ballistic  coefficient  ( com¬ 
pared  to  /3  and  /3*  )  which  results.  In  this  connection,  a  noteworthy  fact  may 
easily  be  established: 

For  each  choice  of  core  and  for  each  coating  with  density 
ratio  p7p<0.5  ,  an  initial  value  of  coating-thickness 

may  be  chosen  so  that  the  ballistic  coefficient  rises  steadily 
toward  p*  as  ablation  proceeds. 

Once  again,  it  is  convenient  to  rewrite  Formula  ( 1)  and  introduce  a  new 
weighting-function  g* 
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^  =  g*(p'.  R)^* 


(8) 


where 

g*(p’.R)  =  C(pV(d(R/r)  +  (l-pVp)(r/R)]  (9) 

for  fixed  r  and  p  .  The  assertion  then  follows  from  the  &ct  that  for  each 
value  p7p<0.5  ,  g*  as  a  function  of  r/R  has  an  absolute  minimum  less 

than  one. 

Propf:—  g*  has  the  form  cx  +  d/x  where  x  =  r/R  , 
d  =  p'/p  is  the  density  ratio,  c  =  l-d>0  ,  and  cd 

is  constant  for  a  given  choice  of  coating.  The  identity 

[cx  +  d/x'f  =  4cd  +  [cx— d/x]^ 

shows  that  an  absolute  minimum  is  attained  when  the  r.h. 
term  vanishes: 

min  r  cx  +  d/x  1  *  /4c3’  =  2cx* 

0<x<l 

where  the  critical  value  of  x  which  gives  the  minimum  is 

X*  =  .  (  End  of  proof. ) 

It  follows  that  g*  has  the  minimum  value 

min  g*  =  /4d  (1  —  dj  (10) 

where  d  =  p'/p  and  that  this  minimum  is  attained  for 

r/R=vmT^d)-  .  (11) 
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Moreover ,  the  simple  identity 


4d(l-d)  =  1  -  (l-2d)^  (12) 

shows  that  the  conditions 

ming*  <  1  with  r/R  <  1  (13) 

+ 

are  compatible  only  for  p'/p  =  d  <  1/2 

An  initial  choice  of  r/R  thus  dictated  by  the  density-ratio  d  =  p'/p 
minimizes  g*  ,  so  different  values  of  r/R  encountered  as  R->r  cause  an 
increase  in  ballistic  coefficient  along  with  g*-»  1  . 

Of  course,  this  argument  shows  that  larger  values  of  r/R<  1  --which 
is  to  say ,  smaller  values  of  R  relative  to  a  fixed  r  --than  that  suggested  by 
Equation  ( 11 )  may  be  chosen  with  the  same  result ,  an  increasing  ballistic  coeffi¬ 
cient.  Equation  ( 11 )  merely  gives  the  critical  value  below  which  a  period  of  initially 
decreasing  P  must  be  tolerated. 

The  corresponding  critical  valiie  of  coating-thickness  is  easily  found  from 
Equation  ( 11) ; 


R  —  r  =  r/'p/p'  —  f  —  r 


(14) 


t  As  a  function  of  two  variables ,  g'*  has  a  saddle -point  at  (p'/p)  =  1/2  ,  (r/R)  =  1. 
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Several  plots  of  g*  are  shown  in  Fig,  2  for  values  of  p'/ p  in  the 
range  0. 01  to  0. 5  .  These  all  lie  above  a  straight  line  from  the  origin  to 
the  point  (  1>  1)  corresponding  to  the  degenerate  case  g*(0,  R )  =  r/R  where 

p’  =  0  . 

The  rise  of  g*  from  its  minimum  is  seen  to  be  almost  linear  for  low 
values  of  p'/p  and  may  be  estimated  for  such  cases  of  interest  by  evaluating 
the  derivative  at  r/R  =  1  : 


£g*(p’.x)i  =  l-2(pVp)  .  (15) 

"x  =  1 

(  This  expression  occurred  once  before  in  Equation  (11)  and  thus  relates  the  rate 
of  growth  to  the  depth  of  the  minimum. )  For  low  values  of  p’/p  ,  /?  thus 

approaches  /3*  at  a  rate  proportional  to  /?*  itself; 


d(f/R)~ 


ri-2(p7p)li8* 


The  critical  values  of  r/R  and  the  corresponding  minima  of  g*  for 
different  density -ratios  p'/p  are  plotted  in  Fig.  3. 
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SUMMARY  AND  CONCLUSIONS 


1.  For  the  usual  case  where  the  coating  has  lower  density  than  the  material 
being  protected,  a  coated  cylinder  of  given  outside  dimensions  has  a  lower  ballistic 
coefficient  than  an  uncoated  cylinder  having  these  dimensions.  The  difference  in 
ballistic  coefficient  is  directly  proportional  to  the  difference  in  densities  and  to  the 
coating-thickness  when  the  latter  is  small. 

2.  The  ballistic  coefficient  of  a  coated  cylinder  can  vary  with  coating  thickness 
in  different  ways  relative  to  the  ballistic  coefficient  of  the  (smaller)  bare  cylinder 
being  protected: 

(a  )  A  coating  whose  density  exceeds  one -half  that  of  the  core  gives  a 
higher  ballistic  coefficient  than  that  of  the  core  alone  for  all  (positive) 
coating -thicknesses . 

(b)  A  coating  whose  density  is  less  than  one -half  that  of  the  core  gives 
a  ballistic  coefficient  which  may  be  either  higher  or  lower  than  that  of 
the  core  and  either  increasing  or  partly  decreasing,  depending  on  the 
range  of  values  of  coating-thickness  considered. 

3.  When  the  density  of  the  coating  is  less  than  one -half  that  of  the  core  ,  initial 
values  of  coating-thickness  can  be  chosen  so  that  the  ballistic  coefficient  rises  steadily 
toward  that  of  the  bare  core  as  ablation  proceeds.  The  density-ratio  determines  a 
maximum  allowable  coating-thickness  for  this  rise  to  occur;  if  a  thicker  coating  is 
used ,  the  ballistic  coefficient  first  decreases  to  a  minimum  which  is  less  than  the 
ballistic  coefficient  of  the  core.  This  minimum  is  also  determined  by  the  density-ratio 
and  corresponds  to  the  critical  value  of  coating-thickness. 
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APPENDIX 


C^A  FOR  A  CYLINDRICAL  SURFACE 

The  central  concept  of  Newtonian  impact -theory  is  the  assumption  that  fluid 
particles  hitting  an  exposed  surface  lose  entirely  their  component  of  momentiun  normal 
to  the  surface  while  the  tangetial  component  remains  unchanged. 

We  denote  by  n  the  inner  normal  to  the  surface  of  the  cylinder  with  the  coor¬ 
dinate  system  shown  below; 


The  free-stream  velocity  v  has  the  direction  of 

00 

U  =  [cos  Q!  ,  0  ,  sinof  ]  ; 

n  =  [0  ,  -sincp  ,  -coscpj  . 

and  the  component  of  normal  to  the  surface  is  therefore  of  magnitude 


V  =  V  .n  =  (-sina  cosmlv 
n  oo  '  oo 
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The  volume  of  fluid  striking  the  element  of  area 


dS  =  Rdcpdx 

in  unit  time  (where  the  cylinder  has  radius  R  )  is  simply  v^dS  ,  so  if  the  fluid 
has  density  ,  the  change  in  the  normal  component  of  momentum  in  unit  time 
corresponding  to  the  mass  -change  dm  is 

2 

V  dm  =  V  (p  V  dS)  =  p  (v  .n)  dS 
n  n  ®  n  ' 

The  resulting  pressure -force  on  the  surface  element  dS  ,  according  to  the 
Newtonian  approximation  is  ,  therefore , 

dff  =  (p  V  sin  Q!  cos  cpRdcpdx)n 

OO  00 

The  component  in  the  direction  of  the  z-axis  is  (cos  cp)dF  ;  symmetry  with  respect 
to  cp  over  the  surface  cancels  out  the  y -components ,  and  the  x -component  of 
is  zero;  integration  over  the  surface  exposed  to  the  flow  thus  leads  to  a  resultant 
force  acting  normally  to  the  cylinder  axis: 

cos^  cp  dcp  dx 


2  2 

F-  =  p  V  R  sin 
N  •  00 


L  37r/2 

;  j 

0  7r/2 


where  R  is  the  radius  ,  and  L  is  the  length. 

I 
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The  integrand  Is  a  periodic ,  even  function,  and  since 


3  2 

cos  cp  =  cos  cp  —  sin  cp  cos  cp  , 


37r/2  ir/2 

^  COSTCO  dcp  =  \  cp  dcp  =  -4/3  • 

ff/2  0 


Thus  ,  in  terms  of  the  dynamic  pressure 


1  2 

q  =  -o  p  V 
^08  2  ®  ® 


Newtonian  impact  theory  predicts  a  resultant  force  (neglecting  end -effects)  represented 
by: 

"n  ■ • 

The  drag  force  is  the  component  acting  in  the  direction  of  the  relative  free -stream 
velocity  v^  : 


Hence,  the  Newtonian  expression  for  the  drag-area  of  a  cylindrical  surface --neglecting 
end-effects --is  given  by  the  formula 


Cj^A  =  I  RLsin^O! 
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A  solid  cylinder  of  density  p  ,  radius  R,  and  length  L,  has  weight 


W  =  irR\p  . 

Its  ballistic  coefficient  accordingly  is 

W/CjjA  =  1.2Rp/sin^a 

where 

3ir/  8~  1. 2  . 


12 


TABLE  1 


APPROXIMATE  DENSITIES 

abs/ft^ ) 


Aluminum 

170 

Nylon 

70 

Beryllium 

115 

Polyethylene 

55 

Chromium 

410 

Sealing  Wax 

no 

Columbium 

525 

Tantalum 

1035 

Copper 

555 

Teflon 

130 

Gold 

1200 

Tungsten 

1200 

Graphite 

140 

Uranium 

1170 
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